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ABSTRACT 

A signal from decaying dark matter (DM) can be unambiguously distinguished from spec- 
tral features of astrophysical or instrumental origin by stu dying its spatial distribution. We 
demonstrate this approach by examining the recent claim of lLoewenstein & Kusenkol d2009l) 
regarding the possible DM origin of the 2.5 keV line in Chandra observations of the Milky 
Way satellite known as Willman 1 . Our conservative strategy is to adopt, among reasonable 
mass estimates derived here and in the literature, a relatively large dark mass for Willman 1 
and relatively small dark masses for the comparison objects. In light of the large uncertainty in 
the actual dark matter content of Willman 1, this strategy provides minimum exclusion limits 
on the DM origin of the reported signal. We analyze archival observations by XMM-Newton of 
M3 1 and Fornax dwarf spheroidal galaxy (dSph) and Chandra observations of Sculptor dSph. 
By performing a conservative analysis of X-r ay spectra, we show the absence of a DM decay 
line with parameters consistent with those of lLoewenstein & Kusenkol (|2009t). For M31, the 
observations of the regions between 10 and 20 kpc from the center, where the uncertainties in 
the DM distribution are minimal, make a strong exclusion at the level above 10a. The Fornax 
dSph provides a ~ 3.3ct exclusion instead of a predicted 4a detection, and the Sculptor dSph 
provides a 3a exclusion instead of a predicted 2.5a detection. The observations of the central 
region of M3 1(1-3 kpc off-center) are inconsistent with having a DM decay line at more than 
20cr if one takes the most conservative among the best physically motivated models. The min- 
imal estimate for the amount of DM in the central 40 kpc of M3 1 is provided by the model of 
ICorbelli et ail (1201 Ol) . assuming the stellar disk's mass to light ratio ~ 8 and almost constant 
DM density within a core of 28 kpc. Even in this case one gets an exclusion at 5.7a from 
central region of M31, whereas modeling all processed data from M31 and Fornax produces 
more than 14a exclusion. Therefore, despite possible sys tematic uncertainties, we exclude the 
possibility that the spectral feature at ~ 2.5 keV found in lLoewenstein & Kusenkol ( l2009t) is a 
DM decay line. We conclude, however, that the search for DM decay line, although demand- 
ing prolonged (up to 1 Msec) observations of well-studied dSphs, M31 outskirts and other 
similar objects, is rather promising, as the nature of a possible signal can be checked. An 
(expected) non-observation of a DM decay signal in the planned observations of Willman 1 
should not discourage further dedicated observations. 



1 INTRODUCTION 



Modern astrophysical and cosmological data strongly indicate that 
a significant amount of matter in the Universe exists in the form of 
dark matter. The nature of dark matter (DM) remains completely 
unknown and its existence presents one of the major challenges to 
modern physics. It is commonly believed that the dark matter is 
made of particlesQ However, the Standard Model of elementary 
particles fails to provide a DM candidate. Therefore the hypothesis 
of a particle nature of DM implies extension of the Standard Model. 



1 Although more exot ic possibilities (such as, for example, primordial 
black holes. lCarrl 120051) 1 exist. 



Weakly interacting massive particles (WIMPs) are probably 
the most popular class of DM candidates. T hey appear in many 
extensions of the Standard Model (see e.g. iBertone et al.l 1 20051 : 
iHoopej [20091) . These stable particles interact with the SM sec- 
tor with roughly electroweak strength jLee & Weinbergll 19771) . To 
provide a correct dark matter abundance, WIMPs should have 
mass in the GeV range. A significant experimental effort is de- 
voted to the detection of the interaction of WIMPs in the Galaxy's 
DM hal o with laborato ry nucleons (direct detection experiments, 
see e.g. iBaudisI | |2007|) ) or to finding their an nihilation s i gnal i n 
space (indirect d e tection experiments, se e e.g. ICarr et alj HoOrj) ; 
iBergstrorrj fe008l) ; lHooper & Baltj l l2008l) ). 

Another large class of DM candidates are superweakly inter- 
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acting particles (so called super-WIMPs), i.e. particles, whose in- 
teraction strength with the SM particles is much more feeble than 
the weak one. Super-WIMPs appear in many extensions of the 
Standard Model: extensions of the SM by right - handed neutrinos 



I by 

( Dodelson & Widrowl [l994l : lAsaka et alj 120051: lLattanzi & Valid 
20071). supersymmetric t heories i Takavama & Yamaguchil l200d : 



Buchmuller et al.1 120071: iFeng et al 

extr a dimensions i Feng et alj l2003h and string-motivated mod 



2004, 2003), models with 



els dConlon & Ouevedq 1 20071) . The feeble interaction strength 
ma kes the laboratory detection of sup er-WIMPs challenging (see 
e.g. iBezrukov & Shaposhnikovl d2007h ). On the other hand, many 
super-WIMP particles possess a 2-body radiative decay channel: 
DM — > 7 + ^,7 + 7, producing a photon with energy E^ — 
Mdju/2. One can therefore search for the presence of such a 
monoch romatic line with energy in the spectra of DM -containing 
objects jAbazaiian et alJl200ll : [Dolgov & Hansenll2002l) . Although 
the lifetime of any realistic deca ying DM is much longer than th e 
lifetime of the Universe (see e.g. lBovarskv & Ruchavskivl d2008h ). 
the huge amount of potentially decaying DM particles in a typical 
halo means a potentially detectable decay signal in the spectra of 
DM-dominated objects. Searching for such a line provides a major 
way of detection of super-WIMP DM particles. 

The strategy of the search for decaying dark matter signal 
drastically differs from its annihilation counterpart. Indeed, the de- 
cay signal is proportional to the DM column density - integral of a 
DM distribution p DM along the line of sight S — J ; pDAi(r)dr 
as opposed to J ; Pr>M( r )^ r m the case for annihilating DM. 
The DM column density va r ies slow ly with the mass of DM ob- 
jects teovarskv et~ai]|2006bl, 12009 ah and as a result a vast variety 
of astrophysical objects of diffe rent nature would produce a com- 
parable decay signal llBovarskv et al]|2006bl : iBertone et alj 120071 
iBovarskv et alj|2009al) . Therefore (a) without sacrificing the ex- 
pected signal one has a freedom of choosing observational targets, 
avoiding complicated astrophysical backgrounds; (b) if a candidate 
line is found, its surface brightness profile may be measured (as it 
does not decay quickly away from the centers of the objects) and 
distinguished from astrophysical lines that usually decay in out- 
skirts of galaxies and clusters. Moreover, any tentative detection in 
one object would imply a signal of certain (comparable) signal-to- 
noise ratio from a number of other objects. This can be checked and 
the signal can either be unambiguously confirmed to be that of DM 
decay origin or ruled out. This allows to distinguish the decaying 
DM line from any possible astrophysical background and therefore 
makes astrophysical search for the decaying DM another type of a 
direct detection experiment. 

In this paper we demonstrate the power of this approach. 
We check the recent conjec ture by lLoewenstein & Kusenkol 
dLoewenstein & Kusenkoll2009l LK09 in what follows), who re- 
ported that a spectral feature at 2.51 ± 0.07 keV with the flux 
(3.53 ± 1.95) x IO -6 photons cm" 2 s _1 (all errors at 68% 
confidence level) in t he spectrum of the Milky Way satellite 
known as Willman 1 dWillman et alj 120051) may be interpreted 
as a DM decay line. According to LK09, the line with such 
parameters is marginally consistent th e restrictions on decayin g 
dark matter from some objects (see e.g. lLoewenstein et alj d2009T) ; 
lRiemer-S0rensen & Hanser] d2009h or the restrictio ns from the 7 ks 
observa tion of Ursa Minor dwarf spheroid al in IBovarskv et al] 
d2007bh). The results of other works (e.g. IWatson et alj 1 20061 : 
lAbazaiian et alj [2007: Boy arsky et alj[2008al) are inconsistent at 
3<r level with having 100% of decaying DM with the best-fit pa- 
rameters of LK09. However, it is hard to exclude completely new 
unknown systematic uncertainties in the DM mass estimates and 



therefore in expected signal, from any given object. To this end, in 
this work we explicitly check for a presence of a line with param- 
eters, specified above, by using archival observations by XMM- 
Newton and Chandra of several objects, where comparable or 
stronger signal was expected. We compare the signals from several 
objects and show that the spatial (angular) behavior of the spectral 
feature is inconsistent with its DM origin so strongly, that system- 
atic uncertainties cannot affect this conclusion. We also discuss a 
possible origin of the spectral feature of LK09. 



2 CHOICE OF OBSERVATIONAL TARGETS 

The particle flux of the decaying dark matter (in photons cm -2 x s)) 
is given by 
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(1) 



where m DM is the mass of the dark matter particle, V is the decay 
rate to photons (depending on the particular model of decaying dark 
matter), is the dark matter mass within the instrument's field- 
of-view (FoV) flf ov and Dl is the luminous distance towards the 
object^ 

Willman 1 signal. In LK09 the signal was extracted from 
the central 5' of the Willman 1 observation. Assuming the 
DM origin of this signal, the observed flux (3.53 ± 1.95) x 
10 -6 photons cm^s -1 corresponds to Fwi = (4.50 ± 2.5) x 
10~ 8 photons cm _2 s arcmiii -2 flux per unit solid angle. The 
mass within the FoV of the observation of Willman 1 was esti- 
mate d in LK09 to be M§ \ = 2 x 10 6 M© , using the best-fit values 
from lstrigari et alj J2008I) . Taking the luminous distance to the ob- 
ject to be Dl = 38 kpc, the average DM column density of Will- 
man 1 is <Swi — 208.5 Mq pc -2 (see however discussion in the 
Section [5] below). To check the presence of a DM decay line, we 
should look for targets with comparable signal. 

The signal-to-noise ratio for a weak line observed against 
a featureless continuum is given by (see e.g. discussion 
in lBovarskv et alj d2007al) ): 



(S/N) oc Ss/t^AeffSliwAE 



(2) 



where 5 is the average DM column density within an instrument's 
field-of-view Qf ov , t C x P is the exposure time, Aes is an effective area 
of the detector and AE is the spectral resolution (notice that the 
line in question has a much smaller intrinsic width than the spectral 
resolution of Chandra and XMM-Newton). For the purpose of our 
estimate we consider the effective areas of XMM-Newton MOS 
cameras and Chandra ACIS-I (as well as their spectral resolution) 
to be approximately equal, and effective area of PN camera is 2 
times bigger[f| 

Galactic contribution. The contribution of the Milky Way's 
DM halo along the line of sight should also be taken into account, 
when estimating the DM decay signalQ Using a pseudo-isothermal 
profile 



Piso(r) 



Pc 



1 + r 2 /r' 2 



(3) 



2 Eq. {T} is valid for all objects whose size is much smaller than Dl ■ 

3 See e.g. http : / / xmm . esa . int / external /xmm_user_support /document at io 
and jhttp : / / cxc . harvard . edu/ ccr/ proceedings / 7_proc /pre sent at ions /< 

4 Notice, that both in LK09 and in this work only instrumental (particle) 
background has been subtracted from the observed diffuse spectra. 
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Figure 1. Comparison of DM column densities of the Milky Way (red 
solid line) with the density profiles of M31, Fornax and Sculptor dSphs, 
discussed in the text. The DM column density esti mate for Willman 1 i s 
based on the "optimistic" DM density profile from IStrigari etal.1 (2008), 
used in LK09 (see Section fSTI for discussion). The angle tj> marks the di- 
rection off Galactic center (for an object with galactic coordinates (I, b) 
cos <j> = cos I cos b). 



with parameters, adopted in iBovarskv et al] J2006bll2007bh (r c = 
4 kpc, p c = 33.5 x 1O 6 M / kpc 3 , r© = 8 kpc), the cor- 
responding Milky Way column density in the direction of Will- 
man 1 (1 20.7° off Galact i c cente r) is 73.9M Q pc~ 2 . As demon- 
strated in IBovarskv etaD d2008bl) . the value of DM column den- 
sity computed with this DM density profile coi ncides with the 
be st-fit Navarro-Frenk - Wh ite dNavarro et alj 19971 . NFW) profiles 
of iKlypin et all d2002l) and lBattaglia et all <2005h within few % for 
the off-center angle <f> > 90° . Other dark matter distributions pro- 
duce systematically larger values of dark matter col umn density. 
For ex ample, the best-fit pseudo-isothermal profile of Ke rins et al.1 
< l200ll) would corresponds to an additional Milky Way contribution 
~ 119M0 pc -2 in the direction of Willman 1. The comparison 
between the DM column densities of several objects with that of 
Milky Way in their directions is shown in Fig.[T] 



3 XMM-NEWTON DATA ANALYSIS 

For each XMM-Newton observation observation, we extracted 
both EPIC MOS and PN spectraQ We select events with PATTERN 
<= 12 for MOS and PATTERN == for PN camera and use 
FLAG == 0. We used SAS task edetect.chain to exclude 
point sources, detected with the likelihood values above 10 (about 
4a). For EPIC PN camera, bright strips with out-of-time (OOT) 
events were also filtered out, following the standard procedure[f| 
We imposed stringent flare screening criteria. To identify "good 
time intervals" we performed a light-curve cleaning, by ana- 
lyzing temporal variabili ty in the hard ( E > 10 keV) X-ray 



band (as described e.g. in Read & Ponman 2003 ; Nev alainen et al.1 



l2005t ICarter & Readll2007l) . The resulting light curves were in 



by more than 2<r). The extracted spectra were checked for the pres- 
ence of residual so ft proton contaminants by comparing count rates 
in and out of FoV jPe Luca & Molendill2004l : lLeccardi & Molendil 
l2008h FI In all spectra that we extracted the ratio of count rates in 
and out of FoV did not e xceed 1.13 which indicates thorough soft 
proton flare cleaning (c.f. lLeccardi & Molendill2008l). 

Following the procedure, described in ICarter & Read! d2007t) 
we extracted the filter-wheel closed background from the same 
region as signal (in detector coordinates). The background was fur- 
ther renormalized based on the E > 10 keV band count rates of the 
extracted spectra. The resulting source and backgro und spectra are 
grouped by at least 50 counts per bin using FTOOL dlrbv. B J200^) 
command grppha. 

The extracted spectra with subtracted instrumental back- 
ground were fitted by the powerlaw model of XSPEC in the en- 
ergy range 2.1-8.0 keV data for MOS and 2.1-7.2 keV for the 
PN camera. This choice of the baseline model is justified by the 
fact that for all considered objects we expect no intrinsic X-ray 
emission above 2 keV and therefore the observed signal is dom- 
inated by the extragalactic diffuse X-ray background (XRB) with 
the powerlaw slope F = 1.41 ± 0.06 and the normalization (9.8 ± 
1.0) x 10~ 7 photo ns cm~ 2 s~*keV~*arcmin~ 2 at 1 keV (90% 
CL errors) (see e.g. ITumb et alj|2002l: |Pe Luc~& Molendi j [20041: 
Nevalainen et al.l2 005; Hickox & Markevitch 2006: Ca rter & Read 



20071 : iMoretti et alj|2009l) . Galactic contribution (both in absorp- 
tion and emission) is negligible at these energies and we do not 
take it into account in what follows. For all considered observa- 
tions powerlaw model provided a very good fit, fully consistent 
with the above measurements of XRB. The powerlaw index was 
fixed to be the same for all cameras (MOS1, MOS2, PN), observ- 
ing the same spatial region, while the normalization was allowed to 
vary independently to account for the possible off-a xis calibration 
uncertainties between th e different cameras, see e.g. iMateos et all 
l2009l : ICarter et alJlioifl The best-fit values of the powerlaw in- 
dex and normalizations are presented in Table Q] 

To check for the presence of the LK09 line in the spectra, we 
added a narrow (AE/E ~ 10~ 3 ) Gaussian line (XSPEC model 
gaussian) to the baseline powerlaw model. Assuming the DM 
origin of this line, we fix its normalization at the level Fwi , derived 
in LK09, multiplied by the ratio of total (including Milky Way con- 
tribution) DM column densities and instrument's field-of-view. We 
vary the position of the line in the interval corresponding to the 3cr 
interval of LK09: 2.30 - 2.72 keV. Eq. (O allows to estimate the 
expected significance of the line. 

When searching for weak lines, one should also take into ac- 
count uncertainties arising from the inaccuracies in the calibra- 
tion of the detector response and gain. This can lead to system- 
atic residuals caused by calibration inaccuracies, at the level ~ 5% 
of the model flux (s ome of them having edge -like or even line- 
like sh apes) (see e.g. | Kirsch et al J 1 20021 . 1 20041) as well as discus- 
sion in IBovarskv et all J2008ch for similar uncertainties in Chan- 
dra. To account for these uncertainties we perform the above pro- 
cedure with the 5% of the model flux added as a systematic error 
(using XSPEC command systematic). 



spected "by eye" and additional screening of soft proton flares 
was p erformed by employing a method similar to that desc ribed 
in e.g. lLeccardi & Molendil d2008l) : lKuntz & Snowdenl d2008l) (i.e. 
by rejecting periods when the count rate deviated from a constant 



5 We use SAS version 9 . . and XSPEC 12.6.0. 

B As described for example at http : / /xmm . esac . esa . int /sas /current^ ctoailAteiatEbttitep^ 



7 We use the script h ttp : //xmm2 .esac .esa. int /external /xmm_sw_cal/backgr 
(20(3) provided by the XMM-Newton EPIC Background working group. 

xmm_ 



_sw_cal /backgroun 
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ObsID 


PL index 


PLnorm, 10~ 4 ph. cm" 2 s" 1 keV -1 ) 
at 1 keV (MOS1 / MOS2 / PN) 


X 2 /dof 


M31on observations: 


0109270101 


1.33 


5.71/5.80/7.88 


1328/1363 


0112570101 




5.37/6.12/6.86 




0112570401 




6.13/6.66/6.37 




M31of f observation: 


0402560301 


1.46 


6.28/6.12/9.97 


895/997 


M31out observations: 


0511380101 


1.53 


2.75 /2.88 /3.89 


875/857 


0109270401 


1.47 


5.02/4.41/5.22 


606/605 


0505760401 


1.27 


1.91/2.04/2.47 


548/528 


0505760501 


1.28 


2.25/ 1.85/3.03 


533/506 


0402561301 


1.57 


2.94/3.10/3.95 


489/515 


0402561401 


1.47 


4.33/3.58/4.47 


751/730 


0402560801 


1.55 


5.22/4.49/5.83 


902/868 


0402561501 


1.44 


2.70 / 4.29 / 6.67 


814/839 


0109270301 


1.39 


4.29/4.58/3.75 


413/436 


Fornax observation: 


0302500101 


1.52 


4.46/4.02/3.78 


938/981 



Table 1. Best-fit values of powerlaw index and normalization for each XMM-Newton observation analyzed in this paper. The systematic uncertainty is 
included (see text). In M3 Ion region, powerlaw index for different observations is chosen to be the same, because they point to the same spatial region. 



4 ANDROMEDA GALAXY 
4.1 Dark matter content of M31 

The Andromeda galaxy is the closest spiral galaxy to the 



ied over the years (see e. 


j. Kerins et al.ll200ll;lKlvDin et alj|2002l; 


Widrow & Dubinskill2005 


;lGeehan et al. 2006; Temoel et alj2007l; 


Chemin et al.l2009l; Corbelli et al.l2010t and references therein) for 



an incomplete list of recent works. The total dynamical mass (out 
to ~ 40 kpc) can be determined from the rotation curve measured 
from Hi kinematics. The major uncertainty in determination of dark 
matter content is then related to a separation of contributions of 
baryonic (stellar bulge and especially, extended stellar disk) and 
dark components to the total mass[j The baryonic mass is often 
obtained from the (deprojected) surface brightness profile (optical 
or infrared), assuming certain mass-to-light ratio for the luminous 
matt er in the bulge and the d isk of a galaxy. 

iBovarskv et alj J2008ah (B08 in what follows) analyzed dark 
matter distributions of M31, existing in the literature at that time 
in order to provide the most conservative estimate of the expected 
DM decay signal. The DM column density 5 in more than 10 



9 Some works also take into acc ount a presence of superma ssive black hole 
in the center of a galaxy (see e.g Widrow & Dubinski 2005) and additiona l 
contribution of a gaseous disk (c.f. lChemin et al.l2009l : lcorbelli et alj2010h . 
Their relative contributions to the mass at distances of interest turn out to 
be negligible and we do not discuss them in what follows. 



models from the works i 


Kerins et al. 1200 ll; 


KlvDinetal.j2002l; 


1 Widrow & Dubinski 2005; 


Geehan et alj2006l: 


Temcel et all200 71) 



turns out to be consistent within a factor of ~ 2 for off-center dis- 
tances greater than ~ 1 kpc (see gray lines in FigO. The most 
conservative estimate o f dark matter column density was provided 
by one of the models of IWidrow & Dubinskil J2005t) (called in that 
work M31B), marked as red solid line on Fig(2] 

Rec ently two new Hi s u rveys of the disk of M31 were per- 
formed dChemin et al] 120091 ; ICorbelli et all l2010h and new data 
on ma ss distribution of M31 became available! 10 l In lChemin et al.1 
d2009h modeling of Hi rotation curve betw een ~ 0.3 kpc and 
38 kpc was perfor med. IChemin et alj d200S ) used R-ba nd photo- 



metric information l lWalterbos & Kennicutt U98l Il988l) to deter- 



mine the relative contribution of the stellar disk and the bulge. 
Based on this information and taking into account foreground 
and internal extinc tion/reddening effects (see section 8.2.2 in 
IChemin et al .1120091 . for details) they determined mass-to-light ra- 



tio, using stellar population synthesis mod els |Bell et all 20031) 
Tdisk — 1.7Tp) (in agreem ent with those of I Widrow & Dubinskil 



(2005); Geehan et al.1 j2006l) ). Aiming to reproduce the data in the 
inner several kpc, they explored different disk-bulge decomposi- 
tions with two values of the bulge's mass-to-light ratio Tb u i ge — 
0.8Tq and Tbuige — 2.2T q (all mass-to-light ratios are in solar 
units). IChemin et al. I d2009h analyzed several dark matter density 



We thank A. Kusenko and M. Loewenst ein in drawing our attention to 
these works ( Kusenko & Loewenstein 2010). 
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Figure 2. Dark matter column density of M31 as a function of off- 
center distance. The gray lines represent models from jKerins et al.ll200ll; 
Kfypin et alj|2002t lKerinsll2004t Iwidrow & Dubinskil l2005r . iGeehan et alj 



200d ;l Tempel et alj|2007l) analyzed in B08. The red solid line repr esents 
the most conservative model M31B of I Widrow & Dubinskil J2005t) . The 
m odel of l Chemin et alj [2009 ) (red dashed line), the maximum disk model 
Kerins J2004|) (blue dashe double-dotted line) and the "minimal" model 



Corb elli et al Uioiol) (red dotted line) (see text) are shown for compari- 



profiles (NFW, Einasto, cored profile). According to their 6 best-fit 
models (3 DM density profiles times two choices of mass-to-light 
ratio) the DM column density in the Andromeda galaxy is higher 
than the one, adapted in B08 (model M31B) everywhere but inside 
the i nner 1 kpc (see red d ashed line on Fig|2j. 

ICorbelli et alj fcoifj) used circular velocity data, inferred from 
the Hi rotation curves, also extending out to ~ 37 kpc. The 
authors exclude from the analysis the inner 8 kpc, noticing the 
presence of structures in the i nner region (such as a bar), associ- 
ated with n on-circular motions. ICorbelli et alj J2010t> also use opti- 
cal data of IWalterbos & Kennicuttl J1988T) and determine an upper 
and low er bounds on the disk mass-to-light rat io u sing the same 
model s jBell & de Jongl 1200 it iBell et alj|2003h as IChemin et all 
d2009h . They obtain possible range of values of the B-band mass- 
to-light ratio 2.5T© < Tdisk < 8T0. These v alues do not take 
into ac count corrections for the internal extinction !"! Corbelli et alj 
d2010h fit the values of Tdisk rather than fixing it to a theoreti- 
cally preferred value. The mass-to-light ratio of the disk of a spi- 
ral galaxy is known to be poorly constrained in such a procedure 
since the contr ibutions of the disk and DM h alo are similar (see e. g. 
discussions inlWidrow & Dubinskil l l2005t) ; IChemin et alj J200II ). 
ICorbelli etaflfeoich analyzed variety of DM models and mass- 
to-light ratios, providing good fit to the data. The mass model 
that maximizes the contribution of the stellar di sk (Tdisk = 8) 
has the core of the Burkert profile jBurkertlll995l) tb = 28 kpc 
(if one imposes additional constraint on the total mass of M31 
within several hundred kpc). The DM column density remains in 
this model remains essentially flat from the distance ~ rs in- 
ward (see FigO. Notice, that the "maximum disk" fitting o f lKerinsJ 



11 The correction for foreground extinct ion was taken int o account 
in ICorbelli etafl feOld) . using the results of lSeigar et all feOOSt) . The uni- 
form disk extin ction was not applied (which explains high values of Tdj s k). 
The authors of ICorbelli et alj 1201 Oh had chosen not to include any uni- 
form extinction corrections due to the presence of a gradient of B-R index 
(E. Corbelli, private communication). 



ObsID 


Cleaned exposure [ks] 


FoV [arcmin 2 ] 




(MOS1/MOS2/PN) 


(MOS1 /MOS2/PN) 


0109270101 


16.8/16.7/15.3 


335.4/336.0/283.6 


0112570101 


39.8/40.0/36.0 


332.9/333.1 /285.9 


0112570401 


29.8/29.9/23.5 


335.6/336.1 /289.5 



Table 2. Cleaned exposures and FoV after the removal of point sources 
and OOT events (calculated using BACKSCAL keyword) of three M31on 
observations. 



d2004l> (blue dashed double-dotted line on Fig|2]l has higher DM 
content. Therefore in this wo rk we will adopt the Burkert DM 
model of ICorbelli et al. as a minimal possible amount of 

dark matter, consistent with the rotation curve data on M31. The 
corresponding DM column density is factor 2-3 lower than the one 
given by the previously adapted model M31B in the inner 10 kpc. 
In what follows we will provide the restrictions for both M31B (as 
the most conservative among physically motivated models of D M 
distribution in M31) and Burkert model of lCorbelli et alj d2010h . 



4.2 M31 central part 

Three observations of the central part of M31 (obsIDs 
0112570401, 0109270101 and 0112570101) were used in 
B08 to search for a dark matter decay signal. After the removal of 
bright point sources, the diffuse spectrum was extracted from a ring 
with inner and outer radii 5' and 13', centered on M310 (see B08 
for details, where this region was referred to as ring5-13). We 
call these three observations collectively M31on in what follows. 
The baseline powerlaw model has the total \ 2 1328 for 1363 
d.o.f. (reduced \ 2 = 0.974). 

Let us estimate the improvement of the signal-to-noise ra- 
tio (O, assuming the DM origin of LK09 spec tral feature. The aver- 
age DM column density in the model M3 IB of lWidrow & Dubinskil 
< l2005t) ) is 5m31oh = 606 M Q pc 2 . The dark matter column den- 
sity from the Milky Way halo in this direction is 74.1 M@ pc -2 . 

The ratio of t cxv x fifov x A c ff of all observations of M31on 
(Table[2} and the LK09 observation is 12.9. Thus one expects the 
following improvement of the S/N ratio: 



(S/N)i 



(S/N), 



606 + 74.1 
208.5 + 73.9 



V12.90 w 8.65 



(4) 



i.e. the ~ 2.5<r signal of LK09 should become a prominent fea- 
ture (formally about ~ 21.6<j above the background) for the 
M31on observations. As described in Section [3] we add to the 
baseline powerlaw spectrum a narrow Gaussian line with the 



normalization fixed at Fm31oh = 



606+74.1 
208.5+73.9 



Fwi 



1.08 x 



10 7 photons cm 2 s 1 arcmin 2 . The quality of fit becomes 
significantly worse (see Fig. [3}- The increase of the total \ 2 due 
to the adding of a line is equal to (23. 2) 2 , (22.9) 2 and (22. 2) 2 for 
lcr, 2(7 and 3cr intervals with respect to the central respectively. 

In addition to the M31on observations, we processed an 
XMM-Newton observation 0402560301, positioned w 22' off- 
center M31 (RA = 00h40m47.64s, DEC = +41dl8m46.3s) - 
M31of f observation, see Table [3] We collected the spectra from 



12 W e adopt the distance to M31 Dj, = 784 kpc jStanek & Garnavicbl 
1 19981) (at this distance 1' corresponds to 0.23 kpc). 
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ObsID 



0402560301 



Cleaned exposure [ks] FoV [arcmin 2 ] 
(MOS 1 / MOS2 / PN) (MOS 1 / MOS2 / PN) 



41.9/42.2/35.2 



405.6/495.4/433.3 



Table 3. Cleaned exposures and FoV (calculated using BACKSCAL key- 
word) of the observation M31of f (obsID 0402560301). The significant 
difference in FoVs between MOS 1 and MOS2 cameras is due to the loss 
CCD6 in MOS1 camera. 

the central 13' circle. Several point sources were manually ex- 
cluded from the source spectra. The rest of data reduction is de- 
scribed in the Section l4~2l The fit by the powerlaw model is 
excellent, the total \ 2 equals to 895 for 997 d.o.f. (the reduced 
X 2 = 0.898). Using the DM estimate based on the model M31B 
we find that the average DM column density for the M31off 
5M3ioff = 388. 6Mq pc -2 plus the Milky Way halo contribution 
74.3Mq pc -2 . The estimate of the line significance is similar to 
the previous Section and gives 



(S/JV), 



388.6 + 74.3 _ 
208.5 + 73.9 1 



4.85 



(5) 



(S/N) W1 

and therefore, under the assumption of DM nature of the fea- 
ture of LK09, one would expect ~ 12.1a detection. After that, 
we add a narrow line with the normalization Fyi3ioS ~ 7.37 x 
10~ 8 photons cm~ 2 s _1 arcmin -2 and perform the procedure, 
described in Section[3] The observation M3 lof f rules out the DM 
decay line origin of the LK09 feature with high significance (see 
Fig. EJ: the increase of \ 2 due to the addition of this line has the 
minimum value A^ 2 — (10. 4) 2 at 2.44 keV, which is within la 
interval of the quoted central value of the LK09 feature. 

We also perform the analysis, adding a line, whose flux is 
determin ed according to the D M density estimates based on the 
model of ICorbelli etal] j2010h . that we consider to be a minimal 
DM model for M31 (as discussed in l4.lt . For this model the col- 
umn density in the central 1-3 kpc decre ases by a factor ~ 3.4 as 
compa red with the value, based on M31B. lKusenko & Loewensteinl 
claimed that in this case the LK09 line becomes consistent 
with the M31on observations. However, our analysis shows that 
the total x 2 increases, when adding the corresponding line to the 
model, by (5.7) 2 for 3cr variation of the position of the line. The 
corresponding increase of total x 2 f° r M31of f region is (2.0) 2 . 
Combining M31on and M31of f observations, one obtains (6.2) 2 
increase of the total x 2 - 

We conclude therefore that despite the uncertainties in DM 
modeling, the analysis of diffuse emission from the central part of 
M31 (1-8 kpc off the center), as measured by XMM-Newton, dis- 
favors the hypothesis that the spectral feature, observed in Will- 
man 1, is due to decaying dark matter. Nevertheless, to strengthen 
this conclusion further we analyzed available XMM-Newton ob- 
servations of M31 in the region 10-20 kpc off-center, where the 
uncertainties in the mass modeling of M31 reduce significantly 
as compared with t he central 5 — 8 kpc (c.f. IChemin et alj|2009 : 
ICorbelli et abfeoiol . see also Section RTTI in particular Fig. [2j. 



4.3 M31 off-center (10-20 kpc) 

We selected 9 observations with MOS cleaned exposure greater 
than 20 ks (Table f4]l. The total exposure of these observations is 
about 300 ks. Based on the statistics of these observations one 
would expect a detection of the signal of LK09 with the signi f- 
icance 12.1a (for M31B) and 11.2a (for ICorbelli et al.l told) ). 



ObsID 


Cleaned exposure [ks] 
(MOS 1 / MOS2 / PN) 


FoV [arcmin 2 ] 
(MOS 1 / MOS2 / PN) 


0511380101 


44.3 / 44.5 / 37.6 


356.2 / 400.4 / 333.3 


0109270401 


38.5/38.4/33.5 


387.8/384.1 /366.7 


0505760401 


26.0/26.0/21.7 


330.5/374.1 /325.0 


0505760501 


23.8/23.8/ 19.8 


344.8/395.2/313.2 


0402561301 


22.7/22.7/20.1 


367.2/419.5/382.6 


0402561401 


39.3/39.3/33.7 


349.9 / 408.3 / 343.6 


0402560801 


42.6 / 42.6 / 36.6 


373.6/435.8/393.7 


0402561501 


38.5/38.5 /34.0 


369.3/421.0/407.5 


0109270301 


24.4 / 24.6 / 22.2 


443.4/439.6/405.0 



Table 4. Cleaned exposures and FoV after the removal of point sources and 
OOT events (calculated using BACKSCAL keyword) of nine M31 observa- 
tions off-set by more than 10 kpc from the center (M3 lout observations). 



The fit to the baseline model is very good, giving total \ 2 equal to 
5931 for 5883 d.o.f. (the reduced % 2 = 1-008). However, adding 
the properly rescaled line significantly reduced the fit quality, in- 
creasing the total x 2 by (12.0) 2 /(10.7) 2 /(10.7) 2 when varying line 
position within la, 2a and 3a intervals (using M31B model) and 
by (11.7) 2 , ( 10.7) 2 and (10.6) 2 for la, 2a and 3a intervals (using 
the model of lCorbelli et~ai]|2010l which gives DM column density 
in M3 lout about 160 - 180M Q pc~ 2 (including Milky Way con- 
tribution), see Fig.|2j. 



4.4 Combined exclusion from M31 

Finally, performing a combined fit to all 13 observations of M31 
(M31, M31of f, M31out) we obtai n the exclusion of more tha n 
26a (using the DM model M31B of IWidrow & Dubinskil j2005h) 
and m ore than 13cr for the minimal DM model of ICorbelli et al .1 
As described in Section [3] in deriving these results we 
allowed the normalization of the baseline powerlaw model to 
vary independently for each camera, observing the same spatial re- 
gion, added additional 5% of the model flux as a systematic uncer- 
tainty and allowed the position of the narrow line to vary within 
2.3 - 2.72 keV interval. 



5 DWARF SPHEROIDAL GALAXIES 

Since Willman 1 is purported to be a dwarf spheroidal (dSph) 
galaxy (but see section [75TTJ, as the next step in examining the hy- 
pothesis of LK09 we compare estimates of the dark mass of Will- 
man 1 to dark masses estimated for other dSph satellites of the 
Milky Way. Specifically we consider two of the (optically) bright- 
est dSphs, Fornax and Sculptor, for which comparable X-ray data 
exist. To characterize the d ark matter ha l os of d Sphs, we adopt the 
general DM halo model o f IWalker et all d2009al) . with density pro- 
file given by 



p(r) 



Pa 



1 



(6) 



where the parameter a controls the sharpness of the transition 
from inner slope lim r _,.o d In(p) /d ln(r) oc —7 to outer slope 
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Figure 3. Spectra of three XMM-Newton observations of M31 (M3 Ion, Sec. 14.21 . The Gaussian line, obtained by proper scaling of the result of LK09 is also 
shown. The top row is for MOS1, the middle row - for MOS2, and the bottom row for PN cameras. The spectra in the left column are for the observation 
0109270101, the middle column - 1 1 2 5 7 1 1 and the right column - for 0112570401. The error bars include 5% of the model flux as an additional 
systematic error. For all these spectra combined together, the \ z increases by at least 22. 2 2 , when adding a narrow Gaussian line in any position between 
2.3 keV and 2.72 keV (see text). 



lira, — >oo d\n(p) / dln(r) oc —3. This mode l includes as specia l 
cases both cored (a = 1, 7 = 0) and NFW jNavarro et"a7]| 19971) 
(a = 7 = 1) profiles. Assuming that dSph stars are spherically dis- 
tributed as massless test particles tracing the DM potential, the DM 
density relates to observables — the projected stellar density and ve- 
locit y dispersion profiles — via the Jeans Equati on (see Equation 
3 of lWalkeretai] ( 12009^ 1. IWalker et alj d2009ij) show that while 
the data for a given dSph do not uniquely specify any of the halo 
parameters individually, the bulk mass enclosed within the opti- 
cal radius is generally well constrained, subject to the validity of 
the assumptions of spherical symmetry and dynamic equilibrium 
and negligible contamination of stellar velocity samples from un- 
resolved binary-orbital motions. 

5.1 Willman 1 

Despite some indicat ions that Willman 1 is a dark-matter d om- 
inated dSph galaxy jMartin et ail 120071 ; IStrigari et al] 120081) . it 
should be noted that the galactic status of Willman 1 remains uncer- 
tain. Claims that Willman 1 is a bona fide galaxy — as opposed to a 



star cluster devoid of dark matter — stem primarily from two consid- 
erations. First, under simple dy namical mod e ls, the stellar velocity 
dispersion (4.3 ± 2.5 km s- MMartin et al.1 ( 120071) 1 of Willman 1 
implies a large mass-to-light ratio (M/Lv ~ 700Mq/Lv,©), in- 
dicative of a dominant dark-matter component like those that char- 
ac terize other dSph gal axies. Second, the initial spectroscopic study 
bv lMartin et al .1 00071) suggests that Willman 1 stars have a metal- 
licity range —2.0 < [Fo/H] < —1.0, significantly broader than 
is observed in typical star clusters. However, both arguments are 
vulnerable to scrutiny. For example, if the measured velocity dis- 
persion of Willman 1 receives a significant contribution from un- 
resolved binary stars and/or tidal heating, then the inferred mass 
may be significantly overestimated. Demonstrations that binary or- 
bital motions contribute negligibly to dSph velocity dispersions 
have thus far been limited to intrinsically "hotter" systems, with 
velocity dispersions o y ~ 10 km s -1 (Olszewsk i et aljfl995t 
lHargreaves et alj[l996t) . Furthermore, given Willman l's low lu- 
minosity, kinematic samples must inc lude faint stars clo se to the 
main-sequence turnoff (see Fig ure 8 of lMartin et all ( 120071) ). These 
stars are physically smaller than the bright red giants observed in 
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Figure 4. The spectra of off-center M31 observation 0402560301 (M31of f, Sec 14.31 . The decaying dark matter signal, obtained by proper scaling of the 
result of LK09 is also shown. The error bars include 5% of the model flux as an additional systematic error. Fitting these spectra together excludes the properly 
scaled line of LK09, at the level of at least 10.4cr (see text). 



brighter dSphs, and thus admit tighter, faster binary orbits. There- 
fore, the degree to which binary motions may inflate the esti- 
mated mass of Willman 1 remai ns unclear. Finally, follow-up spec- 
troscopy by ISiegel et"al] d20Q8|) suggests t hat th e relatively high- 
metallicity stars observed bv lMartin et al.l d2007() are in fact fore- 
ground contaminants contributed by the Milky Way along the line 
of sight to Willman 1. Removal of these stars fro m the sample 
would erase the large metallicity spread reported by iMartin et al .1 
d2008l) . and the remaining stars would have a narrow metallicity 
dist ribution consistent with that of a metal-poor globular clus- 
ter dSiegel et alj|200^) . Nevertheless, we proceed under the as- 
sumption that Willman 1 is indeed a dSph galaxy for which the 
measured stellar velocity dispersion provides a clean estimate of 
the dark matter content. This assumption gives conservative esti- 
mates of the expected DM decay signal in other objects (including 
M31), since any overestimate of the DM column density in Will- 
man 1 will result in an underestimate of the relative S/N expected 
in other objects (see Eq. d2}). For Willman 1 we adop t the line-of- 
sight v elocity data for 14 member stars observed by IMartin et al .1 
d2007l) . These stars, which lie a mean distance of 30 pc in projec- 
tion from the center of Willman 1, have a velocity dispersion of 
ay — 4.3 ± 2.5 km s _1 . To specify the stellar surface density, 
we adopt a Plummer profile, I(R) = 7o[l + -R 2 / r haif]~ 2 > with 
parameters rhaif = 25 pc and Jo = 0.5L© pc -2 as measured by 
IMartin etalJd2008h . Of course the single point in the empirical ve- 
locity dispersion "profile" relates no information about the shape of 



the profile, and these data therefore provide only weak constraints 
on the dark matter content of Willman 1. The solid black line in 
the bottom-left panel of Figure [5] displays the median spherically- 
enclosed mass profile for the general halo model of Equation [6] 
with dotted lines enclosing the region corresponding to the central 
68% of posterior p arameter space from the Markov-Chain Monte 
Carlo analysis (see lWalker et al .1 d2009al) for a detailed description 
of the MCMC procedure). The small published kinematic data set 
(14 stars) for Willman 1 allows for only a single bin in the velocity 
dispersion profile (see the upper right panel in Fig.|5]and compare 
it e.g. with Fornax dSph (the left panel)). This resulting constraints 
on the mass profile of Willman 1 are therefore very poor. The al- 
lowed masses span several orders of magnitude at all radii of in- 
terest and are consistent with negligible dark matter content — the 
lower bound (68% confidence interval) on the mass within 200 pc 
(includ ing all of the Willman 1) is as low as lOOM^ Jstrigari et al .1 



1 20081) obtained for Willman 1 M(< lOOpc) = 1.31JJ x 1O 6 M 



with the DM column density for the best-fit values of parame- 
ters being Swi ~ 2 00Mp) pc~ 2 . The discrepancy with results of 
I Strigari et alJ J2008h may arise from two possible sources. First, 
IStrigari et alj d2008T) use an independent, still unpublished kine- 
matic data set (B. Willman et al., in preparation) for Willman 1. 
Although the global velocity dis persion is sta t istical ly equivalent 
to that which we measure from Martin et al. I d2008h sample, the 
data set used bv lStrigari et alj d2008l) is larger, containing 47 Will- 
man 1 members. Second, IStrigari et al] (2008) adopt more strin- 
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Figure 5. DM mass modeling in Fornax and Will man 1 . Top panels d isplay empirical velocity dispersion profiles as calculated bv lWalker et alj feoQ9al) for 
Fornax, Sculptor and from the Willman 1 data of iMartin et alj j2007l) (for Willman 1, the small published sample of 14 stars allows for only a single bin). 
Overplotted are the median profiles obtained from the MCMC analysis described by [Walker et al] fe009a[) , using the general dark matter halo model |6), as 
well as the best-fitting NFW and cored halo models. Bottom panels indicate the corresponding spherically-enclosed mass profiles. In the bottom panels, dashed 
lines enclose the central 68% of accepted general models (for Willman 1, the lower bound falls outside the plotting window). Vertical dotted lines indicate the 
outer radius of the field of view of the X-ray observations (14' for Fornax, 7.5' for Sculptor and 5' for Willman 1). 



gent priors motivated by cosmological N-body simulations, under 
which p a and r s are s trongly correlated teullock & Johnstorj2007l ; 
iDiemand et alj|2007T) . We have used more general models here in 
order to help illustrate the uncertainties in the mass modeling of 
Willman 1 . However, in the interest of obtaining the most conser- 
vative estimates for the signal we expect to see from other objects 
if the LK09 detection is real, we shall continue to adopt one of 
the largest available estimates of the dark matter column density 
in Willman 1 — namely, the value of .Swi = 208. 5Mq pc -2 that 
results from the modeling o f lStrigari et al] d2008h . 

5.2 Fornax and Sculptor dSphs 

For Fornax and Sculptor, two of the brightest and most well- 
studied dSph satellites of the Milky Way, the availability of large 
kinem atic data sets of ~ 2600 and ~ 1400 members, respec- 
tively dWalker et al]|2009bh . allows us to place relatively tight con- 
straints on the dark-matter column densities. Figure [5] displays the 
empirical velocity dispersion profiles, as well as the fits we obtain 
from the general (Eq. [6jl, NFW and cored ha lo models. Adopting a 
Fornax distance of ~ 138 kpc (Mateo 1998), we integrate the pro- 



Fornax Sculptor 

Pb [A/ P c" 3 ] r s Pb [M Q P c" 3 ] r s 

NFW 0.1 500 0.04 920 

Core 0.57 280 0.48 350 



Table 5. Best-fit NFW (a = 7 = 1 in Eq.(6)) and core profiles (a = 1, 
7 = 0) for Fornax and Sculptor dSphs 

jected density profile out to a radius of 560 pc, which corresponds 
to the 14' region (extraction region of XMM-Newton observation, 
see below) and is similar to the half-light radius of Fornax. For 
Fornax, the best-fitting density profile parameters are given in Ta- 
ble [50 For NFW density profile (V ma x = 17.0km s" 1 ) gives 

13 Note that the fitting para meters p s and r s ar e typically degenerate i n 
the mass modeling of dSphs jStrigari et aT]|2006l ; Fenarrubia et al.ll20O8l) . 
Thus neither parameter is constrained uniquely. The best fits considered 
here generally represent a "family" of models that follow a p s , r s relation 
consistent with the data, but whi ch all tend to have the same bulk mass 
enclosed within the optical radius jPenarrubia et alj2009h . 
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SFom.NFW = 55.2Mg pc~ 2 and the best-fitting cored profile 
(Vmax = 17.4 kms -1 ) gives a slightly lower value 5p or n = 

54.4Mq pc -2 that we adopt for subsequent ana lysis. 

Adopting a Sculptor distance of 79 kpc dMateol [T998I) . we 
integrate the projected density profile over a square of 175 2 pc 2 , 
which corresponds to the 7.6' 2 region (extraction region of Chan- 
dra observations on the ACIS-S3 chip, see below). The best- 
fit parameters for Sculptor are given in Table [5] The best-fitting 
cored profile (V ma x = 19.8 kms ) gives a DM column den- 
sity c?scui,core = 147M0 pc~ 2 and the best-fitting NFW profile 
(Vmax = 19.4km s -1 ), gives <Sscu1,nfw = 140. 3M pc -2 and 
that we adopt for subsequent analysis. Notic e that using the param - 
eters of the "universal DM density profile" of l Walker et alj l l2009al) . 
one would get only slightly higher values. In any case, the DM 
column densities of Fornax and Sculptor are much more tightly 
constrained than that of Willman 1 . The Fornax column density is 
smaller than the expected column density of Milky Way halo in the 
direction of Fornax, which we estimate to be 83. 2Mq pc -2 . The 
estimated DM column density of the Milky Way in the direction of 
Sculptor is 95. 7M© pc -2 (see Section[2]l. 



5.3 XMM-Newton observation of Fornax dSph 

Fornax is one of the most deeply observed dSphs in X-ray wave- 
lengths. We have processed the 100 ks XMM-Newton observation 
(ObsID 0302500101) of the central part of Fornax dSph. The 
data analysis is described in the Section[3] In each camera, we ex- 
tract signal from the 14' circle, centered on the Fornax dSph. The 
cleaned exposure and corresponding fields of view (calculated us- 
ing the BACKSCAL keyword) for all three cameras (after subtrac- 
tion of CCD gaps, OOT strings in PN camera and CCD 6 in MOS 1 
camera) are shown in Table [6] Based on the DM estimates, pre- 
sented in the Section \52\ we expect the following improvement in 
S/N ratio 



(S/N)t 



(S/N)\ 



54.4 + 83.2 
208.5 + 73.9 



V12.12 « 1.70 



(7) 



i.e. we expect a ~ 4.2a signal from Fornax, assuming DM decay 
line origin of the LK09 feature. 

The fit to the baseline model is very good, \ 2 = 955 
for 983 d.o.f. (reduced \ 2 = 0.972). Note that the param- 
eters are consistent with those of extragalactic diffuse X-ray 
background as it should be for the dwarf spheroidal galaxy 



2 keV (see also|Bovarskv et al 


l2006b!.l2007bl:ljeltema & Profumd 


l2008l:lLoewenstein et al.1 20091: 


Riemer-SSrensen & Hansenll2009l). 



After adding a narrow line with the normalization FF orn « 
2.19 x 10~ 8 photons cm _2 s _1 arcmin~ 2 , the \ increases, the 
minimal increase is equal to (5.1) 2 , (4.2) 2 and (3.3cr) 2 for a posi- 
tion of the line within 1<t, 2<r and 3a intervals from LK09 feature, 
respectively (after adding 5% of model flux as a systematic error). 
We see that adding a properly scaled Gaussian line worsens \ 2 by 
at least ~ 3.3(7 (instead of expected improvement of the quality of 
fit by a factor of about 16 if the LK09 feature were the DM decay 
line), see Fig.|6]for details. 

The combination of all XMM-Newton observations used in 
this work (M31on, M31off, M31out and Fornax) provide a 
minimum of 26.7a exclusion at 2.35 keV (which falls into the 3a 
energy range for the LK09 spectral features), after adding a 5% 
systematic error. Within la and 2a energy intervals, the minimal 
increase of x 2 is 29. 2 and 27. 3 2 , correspondingly. Therefore, one 
can exclude the DM origin of the LK09 spectral feature by 26a. 



XMM-Newton 


Cleaned exposure [ks] 


FoV [arcmin 2 ] 


ObsID 


(MOS 1 / MOS2 / PN) 


(MOS 1 / MOS2 / PN) 


0302500101 


53.8/53.9/48.2 


459.1/548.5/424.9 



Table 6. Cleaned exposures of the XMM-Newton observation of Fornax 
dSph 



This exclusion i s produced by using the m ost conservative DM 
model M31B of IWidrow & Dubinskil d2005t) . Using the minimal 
DM model of lCorbelli et alj ( l2010l) . one can exclude the LK09 fea- 
ture at the level more than 14a. 



5.4 Chandra observations of the Sculptor dSph 

5.4.1 The data and data preparation 

For the Sculptor dSph there is one observation of 50ks (Obsld 
9 5 5 5 ) and 2 1 observations of 5 - 6ks (ObsID 4698-4718), see 
Table [7] for details. In all the observations, the centre of Sculptor 
is on the ACIS-S3 chip, and we have restricted the analysis to this 
chip. 

Throughout th e entire analy s is w e used CI AO 4.1 
with CALDB 4.1 fFruscione et~al] l2006h and XSPEC 12.4 
dArnaud et all |2009|) . All the data were observed using the 
VFAINT telemetry mode and thus required reprocessing prior 
to any data analysis in order to take advantage of the improved 
background event rejection based on 5 x 5 pixel islands instead of 
the normal 3x3 pixel. 

We used the wavelet algorithm wavdetect in CI AO to find and 
exclude point sources for each observation. The removed areas are 
very small compared to the field of view, and were consequently 
neglected in the following analysis. Furthermore the observations 
were lightcurve cleaned using \c_clean. The observations 4 7 05 and 
4 712 were flared and left out in the subsequent analysis. The num- 
ber of point sources removed and the exposure time after lightcurve 
cleaning are given in Tab. [7] The variation in the number of point 
sources is mainly due to slightly different observed areas. The ob- 
servations are all centered on the same spot, but there is a variation 
in detector roll angles. Assuming spherical symmetry of Sculptor, 
this has no importance for the further analysis. 



5.4.2 The spectra 

The spectra were extracted from a square region of 7.6' x 7.6' cov- 
ering most of the ACIS-S 3 chip but excluding the edges where the 
calibration is less precise iChandra X-ray Centrej|2009f) . For later 
fitting they where binned to at least 15 counts per bin. 

The spectra of the short observations (4698-4718, except 
4705,4712) were combined into one observation, see Fig. [7] This 
is justified because they fulfill the following criteria: They are ob- 
servations of the same (non-variable) object and thus have similar 
count rates, they are extracted from the same region of the same 
chip, their rmf 's (response matrix functions) are ext remely similar . 
The spectra were added using the FTOOL mathpha dlrbv. B.ll2008l) 
with the uncertainties propagated as if they were pure Gaussian. 
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Figure 6. The spectra of the XMM-Newton observation 0302500101, Sec. 15.21 The error bars include 5% of the model flux as an additional systematic 
eiTor. Fitting these spectra together excludes the properly scaled LK09 line at the level of at least 2a (3.3<r if one restricts the position of the line to the interval 
2.30-2.72 keV) instead of expected improving of the quality of fit by about 4cr (if the line were of the DM origin). 



5.4.3 Background subtraction 

We subtracted only t he p article background as observed with the 
ACIS detector stoweo 14 l Around 2005 there was a change in the 
spectral shape of ACIS-S3 and consequently new particle back- 
ground files were produced dMarkevitch. M.l l 2009T) . We used the 
new particle background files (2005-2009) for both spectra. The 
particle background was normalised in the 10 — 12 keV inter- 
val where any continu ous emission from Sculptor is negligible 
( Mar kevitch et al 1l2003h . For 9555 the effect of changing the nor- 
malization interval to e. g. 4 — 6 keV is less than 3%. For the 
combined observations there is a larger excess at low energies, 
which makes it more sensitive to the choice of normalization in- 
terval and consequently we chose the conservative approach with 
the 10 — 12 keV interval. The larger excess might come from point 
sources, which are unresolved in the short exposures, but resolved 



14 The files can be 

http://cxc.harvard.edu/contrib/maxim/acisbg/ 



downloaded 



from 



and thus removed in the long exposure. The excess is mainly at 
energies below 2 keV and are irrelevant to the present analysis. 

The background subtracted spectra are shown in Fig. [7] with 
the combined observations in black and ObsID 9555 in a lighter 
(green) color. Also shown is the best fit to a power law model and 
the Gaussian expected according to the signal proposed by LK09 
which is further discussed below. 



5.4.4 Constraints 

The total exposure of Sculptor is 162.1 ks which is 1.63 times 
longer than the Willman 1 observation. Sculptor is observed with 
the ACIS-S3 chip, which at 2.5 keV has an effective area that is 
1.23 times larger than that of ACIS-I which Willman 1 is observed 
with. The Sculptor spectra are extracted from a (7.6') 2 square, 
which is 0.74 of the 5' circle for which the Willman 1 spectra were 
extracted. 

(S/iV)sc _ 140.3 + 95.7 



208.5 + 73.9 



VL63 x 1.23 x 0.74 = 1.02 (8) 
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Energy (keV) 

Figure 7. The background subtracted Chandra spectra of Sculptor dSph with the short combined observations in black and ObsID 9555 in a lighter (green) 
color. Also shown the best lit power law model fitted to both observations simultaneously (index fixed to 1.46) and a Gaussian corresponding to the expected 
DM line emission as suggested by LK09. The dip in the residuals is a visual indication of the absence of any line signal at ~ 2.5 keV. 



ObsID 


Obs date 


Exposure" 


Point sources 6 


4698 


2004 Apr 26 


6060 s 


9 


4699 


2004 May 7 


6208 s 


6 


4700 


2004 May 17 


6104 s 


8 


4701 


2004 May 30 


6069 s 


8 


4702 


2004 Jun 12 


5883 s 


7 


4703 


2004 Jun 27 


5880 s 


8 


4704 


2004 Jul 12 


5915 s 


9 


4705 


2004 Jul 24 




6 


4706 


2004 Aug 4 


6075 s 


8 


4707 


2004 Aug 17 


5883 s 


8 


4708 


2004 Aug 31 


5880 s 


8 


4709 


2004 Sep 16 


6091 s 


10 


4710 


2004 Oct 1 


5883 s 


8 


4711 


2004 Oct 1 1 


5727 s 


9 


4712 


2004 Oct 24 




7 


4713 


2004 Nov 5 


6073 s 


11 


4714 


2004 Nov 20 


5649 s 


13 


4715 


2004 Dec 5 


5286 s 


10 


4716 


2004 Dec 19 


6016 s 


8 


4717 


2004 Dec 29 


6073 s 


11 


4718 


2005 Jan 10 


6060 s 


10 


short obs combined 




113169 s 




9555 


2008 Sep 12 


48935 s 


21 


All obs 




162104 s 





Table 7. The analyzed Chandra observations of Sculptor Dwarf. "The ex- 
posure times after lightcurve cleaning using Icclean. 'The number point 
sources found by wavdetect and removed. 



i.e. we expect a significance of 2.5u in case of a DM signal similar 
to the one of LK09. 

Fitting a powerlaw to the two spectra simultaneously over 
the interval 2.1— lOkeV gives an index of 1.54+^ and a normali- 
sation of (1.221°;^) x 1(T 6 photons keV -1 cirT 2 arcmin~ 2 . The 
fit is excellent with x = 1135 for 1078 d.o.f. (reduced x 2 =l-05). 

The ratio [M f D Z/ D l)scu = 1138M / kpc 2 for Sculptor 
which is only 0.605 of the same ratio for Willman 1 (including 
the Milky Way contribution in both cases) giving a flux normaliza- 
tion of 2.03 x 10" 6 photons cm s . We add this Gaussian to 
the power law model varying the position over the energy interval 
2.1-2.8 keV. The quality of the fit worsen (instead of expected im- 
provement) with the lowest increase in \ 2 being 1141 at 2.3 keV. 
At 2.5 keV the increased value is 1 144 which allows us to exclude 
the possibility of a 2.5 keV feature to be a decay line at the level of 
VH44 - 1135 = 3a. 



6 DISCUSSION 

This work provides more than 80 times improvement of statistics 
of observations, as compared to LK09 (factor ~ 5 times larger to- 
tal exposure time and factor ~ 4 improvement in both the effec- 
tive area and the field-of-view). Even under the most conservative 
assumption about the DM column density such an improvement 
should have led to an about 14cr detection of the DM decay line. 
In our analysis no significant lines in the position, predicted by 
LK09 were found. However, in the XMM-Newton observations 
that we processed there are several spectral features in the range 
2 — 3keV (see Table[8]l. This is not surprising, as the XMM-Newton 
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Observations 


Powerlaw index, 


PL normalization 


Line position 


Significance 


Line flux 


Line flux 




best-fit value 


best-fit value, 


[keV] 




best-fit value 


3(7 upper bound 






10 — 6 ph./(cm 2 s arcmin 2 keV) 






ph. /(cm 2 s arcmin 2 ) 


ph. /(cm 2 s arcmin 2 ) 


M31on 


1.32 


1.66/ 1.82/2.39 


2.55 


1.7(7 


7.35 X 10" 9 


2.09 X 10" 8 


M31of f 


1.42 


1.51 / 1.19/2.22 


2.47 


2.3(7 


6.69 x 10~ 9 


3.10 x 10" 8 


M31out 




.../.../... 


2.52 


1.6(7 


3.76 x 10~ 9 


1.16 x 10" 8 


Fornax 


1.37 


0.79/0.59/0.70 


2.19 


2. 60- 


9.63 x 10~ 9 


2.11 x 10~ 8 



Table 8. Parameters of the best-fit powerlaw components (assuming the powerlaw+gauss model) and of the maximally allowed flux in the narrow 
Gaussian line in the interval 2.1 — 2.8 keV. The slight difference between these values and those of Table[T]is due to the presence of gauss ian component 
and the proportionality of powerlaw normalization to field-of-view solid angle (for M31on); for each observation, the corresponding parameters coincide 
within 90% confidence range. The values of best-fit powerlaw parameters for M31out are not shown because they are different for different observations 
(see TablefTJ. 



gold-coated mirrors have Au absorption edge at ~ 2.21 keV 
dKirsch et al J 1 20021 . 1 20041) and therefore their effective areas pos- 
sess several prominent features in the energy range 2 — 3 keV . 
Acco rding to XMM-NewCon calibration report dKirsch et al J 20021 
l2004h . there is about 5% uncertainty at the modeling of the effec- 
tive area of both MOS and PN cameras at these energies. These 
uncertainties in the effective area can lead to artificial spectral fea- 
tures due to the interaction of the satellite with cosmic rays0 In 
particular, the solar protons with energies of few hundreds keV can 
be interpreted as X-ray photons (so called soft proton flares). The 
interaction efficiency of solar protons with the i nstrument is known 
to be totally different from that of real photons dKuntz & Snowdenl 
120081) . In particular, their flux is not aff ected by the Au edge of th e 
XMM-Newton mirrors. According to iKuntz & Snowdenl d2008l) . 
the spectrum of soft proton flares for EPIC cameras is well de- 
scribed by the unfolded (i.e. assuming that the instrument's re- 
sponse is energy-independent) broken powerlaw model with the 
break energy around ~ 3.2 keV, bknpow/b in XSPEC v. 11. 
Therefore, modeling the spectrum that is significantly contami- 
nated with soft proton flares in a standard way (i.e. by using folded 
powerlaw model) will produce artificial residual excess at ener- 
gies of sharp decreases of the instrument's effective area. 

In our dat a analysis we performed both the standard 
flare screening dRead & Ponmarj 120031 ; iNevalainen eTal] 120051 ; 
ICarter & Readl2007 ) that uses the inspection of high-energy (10- 
15 keV) light curve and an additio nal soft proton flare cleaning. 
As found by IPradas & KerrJ d2005h , it is possible to screen out 
the remaining flares, e.g. by the visual inspection of the cleaned 
lightcurve at low and intermediate energies. To provide the ad- 
ditional cleaning after the rejection of proton flares at high ener - 
gies, we followed the pr ocedure of |Pe Luca & Molendj d2004h ; 
IKuntz & Snowdenl d2008h . leaving only the time intervals, where 
the total count rate differs from its mean value by less than 2<j. 

To test the possible instrumental origin of the features, dis- 
cussed in this Se ction, we first performed only a high-energy light 
curve cleaning dRead & Ponmanl 120031 : INevalainen et alj 120051 ; 



ICarter & Read! 12007 ) and determined the maximally allowed flux 
in a narrow line in the energy interval of interest. After that we 
performed an additional soft proton cleaning. This procedure im- 
proved the quality of fit. The significance contours of the resulting 
"line" (for M3 1 out region) are shown on the left panel in Fig. [8] In 
deriving these limits we allowed both the parameters of the back- 
ground, the position of the line and its total normalization to vary 
(the latter from negative values). Finally, we added an unfolded 



15 See e.g. fhttp : / /www, star . le ■ ac ■ uk/ -amr30/BG/BGTable 



broken powerlaw component bknpow/ b to model a contributio n 
of remaining soft proton flares (see e.g. IKuntz & Snowdenl|2008l) . 
The break energy was fixed as 3.2 keV so no degeneracy with the 
parameters of the gaussian is expected. The bknpow/b pow- 
erlaw indices at low and high energies were fixed to be the same 
for different cameras observing the same spatial region, further im- 
proved the quality of fit and made the line detection totally insignif- 
icant (see right Fig.[8]for details). Our analysis clearly suggests the 
instrumental origin of the line-like features at the energy range 2.1- 
2.8 keV. 

As a final note, we should emphasize that the Chandra ACIS 
instrument used in LK09 has the Ir ab sorption edges near 2.11 
and 2.55 keV (see lMirabal & Nietol20l61 for additional discussion). 
They also report negative result of searching for the LK09 feature, 
in agreement with our findings. 



7 CONCLUSION 

In this work we demonstrated that the DM decay line can 
be unambiguously distinguished from spectral features of other 
origin (astrophysical or instrumental) by studying its spa- 
cial distribution. Many DM-containin g objects wou l d pro- 
vide a comparable DM decay signal feovarskv et all l2006bl : 
iBertone et al J 12007 ; Bovarskv et al. 2009ah which makes such a 



study possible (cf. iBovarsky et alj|2009bl : Iden Herder et alj|2009l: 
lRiemer-S0rensen & Ha nsen 2009]). If an interesting candidate line 
is found in an object with an unusually high column density, the 
differences (always within an order of magnitude) in column den- 
sities with other objects can be compensated by longer exposure or 
bigger grasp of observing instruments. 

To illustrate the power of this strategy, we have applied this 
approach to the recent result of LK09 who put forward a hypoth- 
esis about a decaying DM origin of a ~ 2.51 keV spectral fea- 
ture in Chandra observations of the Milky Way satellite known as 
Willman 1 . Although the parameters of decaying DM, correspond- 
ing to such an interpretation, would lie in the regio n, already ex- 
cluded at least 3cr level by several existing works dWatson et al .1 
l2006tlAbazaiian et al.lPi007l ; |Bovarskv et alj|20083) we performed 
a new dedicated analysis of several archival observations of M31, 
Fornax and Sculptor dSphs. 

The exclusions provided by M3 1 are extremely strong (at the 
level of 10 — 26cr, depending on the observed region). The IOct 
bound is obtained using the off-center (10 — 20 kpc) observations 
and a very conservative model of DM distribution with a heavy 
st ellar disc and a very large (~ 28 kpc) DM core (the model 
htmlp f lCorbelli et al.ld2010h . see SectionPTTV Such a model leaves DM 
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Figure 8. Left: A% 2 = 2.3 (dashed-dotted), 4.61 (dashed) and 9.21 (solid) contours (corresponding to 68%, 90% and 99% confidence intervals) for a thin 
gaussian line allowed by the joint fit of M31out spectra used in this work. The resulting line normalizations are shown for an averaged DM column density 
in the M31out region (about 170A/q pc~ 2 , including the Milky Way contribution). The rescaled LK09 line is marked with a cross x; its normalization is 
therefore 4.50 X 10~ 8 X 208 ™ 73 9 = 2.71 X 10- 8 photons cm 2 s arcmin 2 . It is clearly seen that the hypothesis of DM origin for a LK09 line-like 
feature is excluded at extremely high significance, corresponding to more than £S 10<r (see text). Right: The same as in the left figure but with an addition of a 
bknpow/b component to model the contribution of the remaining soft proton contamination (see text). The break energy is fixed as 3.2 keV so no degeneracy 
with the parameters of the gaussian is expected. The bknpow/b powerlaw indices at low and high energies are fixed to be the same for different cameras 
observing the same spatial region. The line significance drops below lrr in this case. 



column density almost unchanged from about ~ 30 kpc inward. 
Combining all XMM-Newton observations of M3 1 at the distances 
1-20 kp c off-center provides more than 13a bound assuming the 
model bv lCorbelli et al and the 26a bound with the model 

M31B of lWidrow & Dubinskil I200I) (with more physically moti- 
vated values of the mass of the stellar disk of M3 1). 

Fornax dSph provides an exclusion between 3.3a an 5a (de- 
pending on how far from its nominal value 2.51 ± 0.07 keV the 
position of the line is allowed to shift) instead of ~ 4a detection 
expected if the spectral feature of LK09 was of the DM origin. The 
bound of 3.3a corresponds to the interval 2.3 — 2.72 keV and the 
5a exclusion corresponds to 2.44 — 2.58 keV interval. 

To summarize, by comparing the strength and position of the 
feature found by Loewenstein & Kusenko with observations of sev- 
eral DM-dominated objects (M31, Fornax and Sculptor dSphs) we 
found that the hypothesis of dark matter origin of LK09 is excluded 
with the combined significance exceeding 14a even under the most 
conservative assumptions. This is possible because of the large in- 
crease of the statistics in our observations as compared with the 
observation used in LK09. 

To change this conclusion, a number of extreme (and not oth- 
erwise motivated) systematic errors would need to be present in 
determination of DM content of M31, Fornax and Sculptor dSph. 
The DM origin of the spectral feature of LK09 would remain con- 
sistent with existing archival data only if the DM amount in M31, 
Fornax and Sculptor are strongly overestimated and/or the mass of 
Willman 1 is grossly underes timated even compared to the best-fit 
model of lStrigari et ail ( l2008h . It should be stressed that in our cal- 
culations of the expected signal in other objects, we have adopted 
for Willman 1 the DM column density that results from the mass 
modeling of Strigari et al. (2008). Our own modeling indicates that 
the mass of Willman 1 is highly uncertain. Even supposing that 
the measured velocity dispersion of Willm an 1 provides a clean 
estimate of the mass, the general models of IWalker et alj d2009ah 
allow the mass of Willman 1 to vary over several orders of magni- 
tude within the region of interest. Furthermore, plausible scenarios 



such as tidal heating and/or contamination from unresolved bina- 
ries in the kinematic sample would cause the Willman 1 mass to 
be over-estimated. We conclude that the column density derived 
from the modeling of lStrigari et al.l d2008l) is approximately an up- 
per limit. In contrast, the large kinematic samples of Fornax and 
Sculptor provide tighter mass constraints under the assumptions 
of our modeling, and given the larger inherent velocity dispersion, 
the modeling assumptions are more secure for Fornax and Sculptor 
than they are for Willman 1 . As a result, we expect that significance 
of our exclusion of dark matter origin of LK09 feature to be very 
conservative. 

Small line-like features (with the significance ~ 2a) are 
present in the XMM-Newton spectrum at energies between 2 and 
3 keV (see Fig. [8] left and Table [8] for details). Their most prob- 
able cause is the interaction of the satellite with soft protons (see 
discussion in Section|6j which can lead to an appearance of artifi- 
cial spectral features due to the uncertainties in the effective area 
at these energies (c.f. lKirsch et al ] |2002l. 120041) . Indeed, the signif- 
icance of the detection of these features drops below la when the 
residual contamination with soft protons is modeled by an unfolded 
broken powerlaw (see previous Section for details). 

The search for decaying DM is a well-motivated task, having 
crucial importance not only for cosmology and astrophysics, but 
also for particle p hysics. X-ray astronomy has a significan t poten- 
tial in this respect jAbazaiiarl2009l ;ld en Herder et alj2009h . Due to 
the small amount of baryons present in dwarf spheroidal galaxies, 
the uncertainties in DM determinations are smaller (as compared 
to e.g. spiral galaxies) and DM modeling provides tight bounds 
on the amount of DM. However, the above-mentioned uncertain- 
ties make Willman 1 (or other ultra-faint dwarfs) to be bad ob- 
servational targets from the point of view of searching for decay- 
ing DM. On the cont rary, the "classical" dS phs provide excellent 
observational targets ( iBovarskv et"ai]|2006bl) . A signal in Fornax, 
Sculptor (or other "classical" dSphs such as Draco or Ursa Minor 
not considered here because of the absence of archival observa- 
tions for them) would therefore have the potential to provide much 
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more information about the nature of DM than would a signal from 
Wil lman 1. To improve sig nificantly the best available bounds (see 
e.g. iBovarskv et al.ll2009bh and to probe the most interesting re- 
gion of parameter space, an extended (300-500 ks) observations 
of the well-studied dSphs are needed. However as we are search- 
ing for the weak diffuse signal, a special care should be taken in 
choosing the instruments and observational times in such a way 
that m inimize possible flare contamination (see e.g. lBovarskv et all 
l2007bh . Drastic improvement in the decaying DM search is pos- 
sible with a new generation of spectrometers, ha ving large field 
of view and energy resolution close to several eV dBovarskv et all 
l2007al : |PiroetaU 2009; den Herder et alJl2009l ; lAbazaiianll2009Tl 
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